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We introduce three deformations, called a-, fl- and 7 -deformation respectively, of a A-body 
probabilistic model, first proposed by Rodriguez et al. (2008), having g-Gaussians as A —>■ 00 
limiting probability distributions. The proposed a- and /^“deformations are asymptotically scale- 
invariant, whereas the 7 -deformation is not. We prove that, for both a- and /3-deformations, 
the resulting deformed triangles still have g-Gaussians as limiting distributions, with a value of q 
independent (dependent) on the deformation parameter in the a-case (/3-case). In contrast, the 7 - 
case, where we have used the celebrated Q-numbers and the Gauss binomial coefficients, yields other 
limiting probability distribution functions, outside the g-Gaussian family. These results suggest that 
scale-invariance might play an important role regarding the robustness of the g-Gaussian family. 


I. INTRODUCTION 


It is well known that the (properly centered and 
rescaled) sum of N independent (or weakly dependent) 
random variables with finite variance approaches a Gaus¬ 
sian distribution in the A —>■ 00 limit [7] . This fundamen¬ 
tal classical result, known as the Central Limit Theorem 
(CLT), is at the basis of the ubiquity of Gaussian distri¬ 
butions in Nature. The classical CLT, however, cannot 
be applied to a set of strongly correlated random vari¬ 
ables. Therefore, in the context of nonextensive statis¬ 
tical mechanics [5H], it has been argued the existence of 
a generalized CLT for random variables correlated in a 
specific way [331IM] , called g-independence. Alternative 
CLTs, based on a different kind of correlations [33] or on 
exchangeability [5], have also been proposed in the liter¬ 
ature. In all these theorems, the (properly centered and 
rescaled) sum of N correlated random variable has, in 
the A —>■ 00 limit, a q-Gaussian distribution Pq{x), 


/ 9-1 ^( 5 - 1 ) 

V - r(|^)"9 

9e (1,3), 

1 p-r" 

9=1, (1) 

/ 1-Q 3 -q r(^r^) 3,2 

V - 2 r(^)®9 

g G (- 00 ,1), 


where g is a real parameter depending on the nature of 
the correlations, and the q-exponential function is de¬ 
fined as follows: 


:= 

Cq 


[1 + (1 - q)x] + 


9 = 1 , 
9 7^ 1, 


( 2 ) 


with [a;]+ := x9(x), 6(x) Heaviside function. In what 
follows, we will use also the inverse function of the g- 
exponential, the q-logarithm, defined as 


logq(a;) 


ln(a;) 

1-9 


9=1, 
9 7 ^ 1, 


X > 0. 


(3) 


Like Gaussian distributions, g-Gaussians also are ubiq¬ 
uitous in Nature. Indeed, analytical, experimental and 
numerical investigations in biology [ 33 ], economics pilSO]. 
high energy physics m, anomalous diffusion processes 
DEI], dynamics of many-body classical Hamiltonian sys¬ 
tems DEHmill], cold atoms [SlllSlIle], dissipative and 
conservative low dimensional maps |25[ I26j . turbulence 
[3] among other^ have shown that g-Gaussian distribu¬ 
tions appear in the probabilistic analysis of many systems 
in which long-range interactions are present, or ergodic- 
ity lacks. These evidences strongly support the existence 
of a generalized CLT involving g-Gaussians. 

To investigate the conditions under which such a gen¬ 
eralized CLT holds, analytically solvable probabilistic 
models yielding g-Gaussian limiting distributions are of 
paramount importance. In particular, in |20j a proba¬ 
bilistic model for A correlated binary random variables 
was introduced, generalizing the celebrated Leibniz tri¬ 
angle [13]. The proposed model preserves the scale- 
invariance property (see Ref. |20] and below for a defi¬ 
nition) of the Leibniz triangle and it can be rigorously 
proved that this model has g-Gaussians with g € [0,1] 


^ For a regularly updated bibliography on nonex¬ 
tensive thermostatistics and related topics, see 

http://tsallls.cat.cbpf.br/biblio.htm 
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as limiting distributions for TV —> oo. Subsequently, it 
was shown [211 [22] that particular d-dimensional scale- 
invariant probabilistic models with d > 1 have, as lim¬ 
iting distributions, Dirichlet distributions for d > 1, 
whereas for d = 1 q-Gaussians were obtained. The 
ultimate relationship between scale-invariance and q- 
Gaussianity, i.e. the appearance of g-Gaussians as proba¬ 
bility distributions for statistical models in the thermody¬ 
namic limit, is not yet completely clear. In Refs. EHlEa 
it was already conjectured that (asymptotic) scale invari¬ 
ance could be possibly a necessary but not sufficient con¬ 
dition, for the emergence of q-Gaussians. Indeed, in (T^j 
two scale-invariant probabilistic models which are differ¬ 
ent from the one analyzed here, are analytically stud¬ 
ied, showing a limiting distribution different from a q- 
Gaussian. 

With the aim of shedding further light on this problem, 
we address here the robustness under small perturbations 
of the general one-parameter family of scale-invariant 
probabilistic models introduced in [20j . In particular, we 
investigate the stability of the g-Gaussian family in the 
space of probability distributions: this property is indeed 
fundamental for the existence of a generalized GLT yield¬ 
ing g-Gaussian distributions in Nature. In this context, 
exactly solvable probabilistic models are essential tools 
for a rigorous study of the type of correlations and the 
properties required for such generalized theorem. 

We consider two new families of asymptotically scale- 
invariant triangles, namely the a-triangles and /3- 
triangles, which generalize the aforementioned family. 
These deformations are based on the introduction of two 
classes of real numbers, the a-numbers and P-numbers 
respectively, in the same spirit of the Q-numbers, typical 
of the Q-deformations of Lie groups and algebras m- 
As we shall see, despite the deformation, the limiting 
distributions remain g-Gaussians, but with a value of g 
which might depend on the perturbation strength. 

To the best of our knowledge, this is the first article ad¬ 
dressing, for specific probabilistic models, the important 
problem of the robustness of g-Gaussians as attractors, 
a fundamental property involved in the existence of a 
generalized GLT. 

The aforementioned deformations may be considered 
as nontrivial ones, since there is no a priori guarantee 
that an arbitrary deformation should preserve the same 
g-Gaussian behavior for large values of N. As a coun¬ 
terexample, we introduce and study here an alternative 
deformation, that we call ^-deformation, based on the 
classical definition of Q-number |13j used in combina¬ 
torics. This deformation does not generically preserve 
g-Gaussian forms for the limiting distributions. Since 
scale-invariance is violated by the y-deformation, in con¬ 
trast with the a- and jd- ones, a possible conjecture 
might emerge on the necessity of (at least asymptotic) 
scale-invariance for g-Gaussianity (see also |52]L 


II. PRELIMINARIES: LEIBNIZ-LIKE 
TRIANGLES AS PROBABILITY MODELS 


Let us consider a system with N identical elements, 
whose states are characterized by binary variables Xi € 
{0,1}, i = 1, Let us introduce also the proba¬ 

bility rjv^n of having a given configuration {a;i,..., xn} 


with The probability PN,n of having any 


configuration such that Xi = n \s obtained taking 


into account the proper degeneracy 

. n. 


PN,n ■= 


rN,, 


(4) 


The set of values rjv,n can be organized in a triangle, 
in such a way that the element r 7 v,n is the nth element 
of the TVth row. We require that the following Leibniz 
triangle rule (or scale-invarianee property) holds: 

XN,n+l T = rj^—i n- (5) 

The Leibniz triangle r^^\ [T^ can be constructed using 
the Leibniz rule and defining g as follows: 


,(i) ._ 

’ N,0 - 


1 

N + 1 


^( 1 ) _ 
' N,n 


(^JN + 1 


( 6 ) 


The A —>■ oo limiting distribution is the uniform distri¬ 
bution. Gonsidering instead rjv,o = j P G (Oil); the 
A — >■ 00 limiting distribution is a Gaussian distribution, 
being in this case pN,n = ~ he., the bi¬ 

nomial distribution. 

The aforementioned Leibniz triangle has been general¬ 
ized in [20]: for G IN, 

(j,) r(2l/)r(A-I- I^) X^M+ 2 (v-l),n+v-l 

^^■0 r(n)r(A + 2u) ^ " rW ,, , ' 

2 (i ^—1 

(7) 

The triangle (|^ is recovered as the u = 1 particular case. 
Remarkably, it has been proved (20] that, for A —>■ oo and 


n = 2,3,.. 


G) 

Pn.u ■= 


we 


A' 


hav^ 


N,i 


iy»i 2 -\/iz — 1 

A 


7}ji(i/) ( 2 :), 


where Pq^i^y){x) is a g-Gaussian with g = qi{v), 

1 


gi(z/) = 1 - 


V- 1 


< 1 : 


( 8 ) 


(9) 


the subindex 1 will become transparent later on. In the 
n ^ 00 limit, rj^g —>■ 2 ~^ and qi{iz) —>■ 1 as expected 
from the GLT. 


^ Observe that in (^)^iv result can be proved 

using the identity ^,^(^)~^ = Jq — x)'^~^dx and the bi¬ 

nomial theorem. 
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III. ASYMPTOTICALLY SCALE-INVARIANT 
DEFORMATIONS OF THE GENERALIZED 
LEIBNIZ TRIANGLE 


In this Section, we introduce two asymptotically scale- 
invariant deformations of the probabilistic model ana¬ 
lyzed in [20] ■ 


A. The a— numbers and the a— triangles 


The basis of our construction is the notion of a- 
number. 



Definition III.l. Given n € IN U 0, and a > 0, a ^ 1, 
an a-number is the real number defined as follows: 


{n}a := (n -1-1)1- 


1 — a 
1 - 

where we have introduced (n > 1) 


nan{a) n > 1, 

0 n = 0 , 

( 10 ) 


,, n-|-l/ 1 —a\ Q_>.i , . 

«„(»):= —( 11 ) 

The previous definition is such that 

{n}i = lim{n}a = n, (12) 

a—)-l 

SO in the a —1 limit we recover the ordinary numbers. 

The following is a generalization of the factorial of a 
natural number. 

Definition III.2. Given n S IN U 0, we shall call a- 
factorial the number defined as 


{0}q! := 1 and {n}a! — ]^{^}a = ^' cLk{o^), n € IN. 

(13) 






The ordinary factorial is recovered in the a —)■ 1 limit, 
{n}i! = limQ_j.i{n}Q,! = n\. We define now an extension 
of the binomial coefficients. 

Definition III.3. Given the nonnegative integers N and 
n < N, the a-binomial coefficient is defined as 


{N}J 


N 

'^Ja {n}J.{N-n}al 

{ 1 if n = 0 or n = V, 


a. Deformation of the Leibniz-like triangles using the 
a-numbers We want now to deform the family of tri¬ 
angles obtained in [201 using the aforementioned a- 
numbers. Let us start introducing the Leibniz-like a- 
triangle as 


,(i) ._ 

' N,n,< 


1 


n = 0,..., N, (15) 


{N + lU{^h 

which is related to the original Leibniz triangle as 


„(i) = ,,(i) „(i) 

' N.n.a r^N.n.Oi N,n'> 


(16) 


with 


■= 

^N,n,cx ■ 


1 


1 


a-N+iict) 

n N—n ak{a) 
fc=l 


, a7V-|-l(Q:) 


for iV = 0, 

for iV > 0 and n = 0 or n = N, 

for TV > 0 and n = 1,...,TV — 1. 

(17) 


In addition, limQ_>i= 1. As an illustration, the 


a = I instance of family (151 is given by 



,(i) 

N,n, 5 


The a-binomial coefficients share with the Pascal co- 
efficients the property { ^ }^ = { }^, Vo. Again, we 

recover the ordinary binomial coefficients in the a —>■ 1 
limit. In the same fashion as the Pascal triangle, the a- 
binomial coefficients can be displayed forming a Pascal 
a-triangle: 


The product of the Leibniz a-triangle by the Pascal 
triangle (which takes into account the degeneracies) does 
not generically yield a set of probabilities, since 
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Nevertheless we can circumvent this difficulty by prop¬ 
erly renormalizing the triangle to get a new one with 
coefficients 


p(i) 

' N,n 


,( 1 ) 

' N,n 


(^) ' 

2^k—0 [ k J^IV,k,a 


and associated probabilities 


p(i) 

PN,n,a — 1 77 , ' ' N^n.a' 


(18) 


(19) 


The normalized Leibniz a-triangie for cx = ^ looks like 



,di) 

N.n 


Following Ref. [20], we can define a two-parametric 


family of triangles from the Leibniz a-triangle (151 for 
a > 0 and z/ S IN: 


'"Nn 


,( 1 ) 

A^+2(l'— l),n+i/—l,a (i/) (ly) 

Cl I “ ^^N,n,a^N,n^ 


( 20 ) 


* 2 {u— 1 ),u— 1 ,q. 

where, for i/ > 1, we hav^ 

MN,n,a ■= 

, . N-n-kv—l , . u—l I . 

a2u-i{oi) j-j- ak[a) jj- ak+v-i(Oij 


^AZ-t- 2 i^ —1 (Q^) ^^ ^^ ^k(^) 

_ a 2 ^-i(a)(a - 1 ) - l) 


He 


,N—n+i^ _ — 


a" 


1 ) 


(N-n + iy)(n + iy) W afc+,.-i(a) , . 

7V + 2Z. ak(a) ' ^ ^ 

As before, normalization is needed to obtain the family 


N,n,a 


_ _ 

'’AT.n.a — AT fN\M ’ 
A^k—0 \ k / N,k,ct 


(23) 


^ To obtain Eq. \22\ we used the following identities for G IN' 
and a, 6 G R"*": 

1 —{1 l — a ^ 

llfc = l - i-af^ + b) - • 

( 21 ) 


n K fc+g+l l+g+fC 

k = l k-\-a 1+a ’ 


with associated probabilities 
Pn’u,^ ■= 


s(‘") 

N,n,Oi 


(24) 


trivially satisfying T,n=oP^,a = 1 - 

Observe that triangle ( |2^ does not generically fulfill 
the scale-invariance condition ([^. Nevertheless, it is 
asymptotically scale-invariant, i.e.. 




lim 

N^OO 


4u) 

'’W-l.n.a 


iV-!-oo A'') 

-^=?7 fixed ' N,n,a ' ' N,n-\-l,(y. 


= 1 . 


(25) 


This follows from the fact that the normalization 
J 2 n=o n a ^ power-law scaling behavior (as we 

will see later by evaluating p large N), and 

from the limit 


,(i) 

' N-l+2{v-l),n+v-l,, 


N —>-oo 


„(i) , „(i) 

' ^+2(1^—l),n+i/—l,a ' ' A^+2(i/—l),n+iv-jC 


> 1 , (26) 


obtained using the limit qn 


N—^OO 


¥ min{l, a}. 


b. A theorem on the robustness of Leibniz-like a- 
triangles In this Section, we prove that the family of 
deformed triangles (23) still possesses a g-Gaussian as 


limiting distribution for N ^ oo, after proper centering 
and rescaling. However, the limiting value of q is different 
from qi{v), although not dependent on a. 


Theorem III.l. The family of triangles (20) with v € 

(27) 


( 1 , -l-oo) and a S \ { 0 } satisfies the property 
N j,A Ar»i 




»(!/) 

'PVn 




where we have introduced the properly centered and 
rescaled variable 




In Eq. (27), is a q-Gaussian with 

for 0 ^ 1 , 


9„(j^) := 1 - 


1 


1 -i 


- Sa,i ll-p3T fora=l. 
See Appendix [A] for the proof. 


In Fig. I we plot some numerical results both for pf^ 
and for tne g-logarithm 


(28) 


(29) 


A) 


wH'Hx) := 




(x) 


PqAAi^) 


= —X 


(30) 


comparing with the theoretical predictions. Observe that 
the following relation between qa{v) and qiiy) holds: 


1 


1 


l-qa{^) l-qi{v) 


-Hi. 


(31) 
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FIG. 1: Probability distributions (24) and Wa\x) := p‘'°('')(o) v = ‘2, N = 10^ and different values of a; 

the derived asymptotic behavior is also depicted as continuous line for a ^ 1. The corresponding values of ga(2) are 

given by Eq. (291. 


c. Entropy Let us now focus on which entropic func¬ 
tional is extensive for the above model. A natural candi¬ 
date is in principle the nonadditive entropy m 


q{u,Oi) 

^qent 


1 - 

^ Z^n=0 


y 

\n) y N,n,aJ 


1 Pent 


(32) 


Using a generalized entropic form [9l [23l [24] , including 
the nonadditive entropy as particular case, in m it has 
been shown that for a wide class of triangles, the number 
of microstates 0, as a function of the system size N, in¬ 
creases according to the law Q{N) = 2^. This leads to a 
scenario in which the only possible value of gent making 
the entropy (32) extensive for a = 1 is gent = 1, which 
corresponds to the Boltzmann-Gibbs entropy. Similar 
arguments, based on the Laplace-de Finetti theorem, 
also yield gent = 1 [II|- 

The same kind of reasoning of [101 HD applies for the 
class of models we address here, and, therefore, also here 
we expect gent = 1- This has been confirmed by a numer¬ 
ical analysis. In Fig. the g-entropy (32) is plotted as 
a function of N for the particular case v = 2 and a = ^. 
In agreement with the above, we find that the value of 
gent which makes the entropy extensive is indeed 

gent — 1- 



versus N for different values of gent- 


shall call j3-numbers the real numbers defined as follows: 

f 0 if n = 0, 

(^1 - -t 1 = n6„(/3) ifn>l,*'^^^ 

where we introduced 


B. The /3—numbers and the /?—triangles 


bn{P) 1 -|- 


1 

n 


1-/3 
1 - /3"' 


(34) 


Let us consider another deformation of the generalized 
Leibniz triangles introduced in Section |TT| 


Note that, so defined, the /3-numbers are related to 
the a-numbers as 


Definition III.4. Given n G INUO, and /3 > 0, /3 ^ 1, we 


[n]p = {n- 1};3 + 1- 


(35) 






















6 


It follows that [n]i := lim^_>i[n ]/3 = n, hence the /3 —t 1 
limit recovers the ordinary numbers. The /3-factorial can 
now be defined as 


[ 0 ]^! := 1 , [n]/3\-.= Y[[k]i3 = n\Y[h{l3), n > 1 . 






(36) 

The factorial number is recovered as [n]i! := 

lim^_>.i[n]^! = n\. We can further define the /3-binomial 
coefficient as 




. /3 [n]i}\[N - n]j3\ 


— 1 -r^N—n bi~ 


'N\ TT-^~ 

‘■nJ llfc=l bk{l3) 


(/ 3 ) 


if n = 0 or n = 
if 1 < n < iV — 1. 


(37) 


where, as expected, lim; 3 _n = (^). 

We can introduce therefore the Leibniz-like /3-triangle 
as 


,(i) 

N,n,P 


1 


1 




n = 0,..., N. 


(38) 


The Leibniz triangle is obtained in the /3 —>■ 1 limit, 
lim/ 3 _n r^^^^ ^ As before, in order to get a set of 

probabilities we have to normalize the triangles, obtain¬ 
ing 


p(i) 

' N,n,l3 


^( 1 ) 

' N,n,P 


z^/c=0 ' ' " » 


(39) 


Jk=0 \ kr N,k,/: 

whose associated probabilities 


j5(i) 

PN,n,l3 


N 


di) 

N,n,l3 


(40) 


satisfy by construction the normalization condition 

p(i) _ I 
l^n=Qt>N,n,P ~ 

We will now generalize the Leibniz /3-triangle by prop¬ 
erly substracting subtriangles of it. In analogy with 
the previous deformation, let us now introduce a two- 
parameter family of triangles 


' N,n,P - 


^( 1 ) 

' A^+2(i/' —l),n+i/' —1,/3 






where iz S IN and 


(.) 62.-i(/ 3) fr'6fc+._i(^) 


bN+2.-l{P) 


n 


hkiP) 


n 

k^l 


(41) 


bk{l3) 


■ After the needed normalization, we obtain the fam- 


' N. 


ily p and the corresponding probabilities 


4’') 

^N,n,P - 


'-M 

PN,n,P - 


^N,n,l3 




E l\ 

k=> 


N 


0 \ kJ' N,k,13 
'’N,n,p- 


(43a) 

(43b) 


It can be proved that the triangle (43a I is asymptotically 
scale-invariant, 


lim 

N—>-oo 


.{u) 

^N-l,n,l3 


jv-i-oo A'') 

-^=?7 fixed N,n,0 ' 


= 1 . 


(44) 


The /3-numbers appear as a variation of a-numbers 
and, moreover, they have the same asymptotic behavior, 
lim„_>oo = 1. Therefore, it could be expected that 

the behavior of a i® same of ^ for A^ —>• 
oo with respect to the parameters of the deformation. 
However this is true only when we consider values of the 
parameters greater than one. Indeed, for /3-triangles, the 
following theorem holds: 


Theorem III.2. The family of triangles (43) with iz € IN 
and /3 > 0 satisfies the property 


N 


(u) p (P 

PN,n,l3 


2\/i^-x(/3) 

where we have introduced the function 

'o 


(45) 


x(/3) := I-Hd/ 3 .i-max |l, ^| = | I for 


for /3 > 1, 

/3 = 1 , 


and the properly centered and rescaled variable 




1-3 /or 0 < /3 < 1, 

(46) 

(47) 


In Eq. (45), is a q-Gaussian with 

I 1^"^ 

qp{i2) = I - 


^z-x(^) -1 


I - 


hi'+i-t 


for 13 > I, 
for /3 = I, 
for 0 < /3 < I. 

(48) 


See Appendix for the proof. 


In Fig. we present some numerical results for ,0- 
triangles for different values of /3. We plot also 


bk+n+ix—l (d) 


(42) 

The /3 —>■ 1 limit of the two-parameters family of triangles 
(41) yields the undeformed family ([^, lini/ 3 -).i = 


w^^\x) := logq^(^) 


PqfiG) (^) 

( 0 ) 


p, 


qpG) 


( 49 ) 


for zz = 2 and different values of /3. 
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FIG. 3: Probability distributions (43bl and := logq^(^) iov i/ = 2, N = 10^ and different values of /?; 




(x) 




the derived asymptotic behavior is also depicted as continuous line for each value of /3: for the sake of clarity not all 
data are plotted. The corresponding values of qp{2) are given by Eq. (48). 


Observe that qa{v) = qp{v) when a > 1 and /? > 1; 
moreover we can write the following relation between the 
limiting value q/ 3 {v) for the deformed triangles (/3 ^ 1) 
and the limiting value qi{v) for the undeformed triangle: 


1-TT “ i-TT + ^ ~ x(/5)- 

l-qi3{v) l-qi{u) 

In particular 

lim -- , , = -- , ' + 1 . 

p^i±l-qp{iy) l-qi{v) 


Eq. (50) can be written, for /3 1, as 


min{/3,1} min{/?, 1} ^ 

l-qpiv) 1-qiiv) 


(50) 


(51) 


(52) 


the previously introduced a-numbers and /3-numbers. 
Eor this reason, in the remainder of the paper we will 
refer to the Q-numbers as the 'y-numhers. 


Note that y-numbers ( |53| ), a-numbers (10) and /3- 
numbers (33) are related as follows: 


n] 


{n 


l}-y + 1 = n 



+ 1 . 


(54) 


In the limit 7 —>■ 1 the 7 -numbers reduce to the ordinary 
numbers, lim.y_j.i|n].y = n. Moreover, a ^-factorial can 
be defined as 


Einally, following the same arguments adopted for the 
a-triangles, it can be easily verified that also in this case 
we have gent = 1 - 


IV. A NON ASYMPTOTICALLY 
SCALE-INVARIANT DEFORMATION 


IOI 7 ! := 1, H 7 ! := ]J[fc]^, linilnlT,! = n\, n> I, 

A A 7—>-1 

(55) 

as well as the standard Gauss binomial coefficients 


Inspired by the Q-calculus [13], we shall consider an 
alternative deformation of the Leibniz triangle based on 
the so called Q-numbers, defined as: 

M 7 — 1 7 G (0,oo) \ {!}, (53) 

1-7 


[VI 7 ! 


|nl7!|iV-nl^! 


7 ,i-j,7 

f 1 


y-fN-n 1 - 7 '“+" 

Ilfc=l 


1 - 7 '“ 


if n = 0 or n = V, 
if 0 < n < 3V, 


(56) 


where we have used the notations 7 and |n].y instead of 
the usual ones in order to avoid confusion with the en- 
tropic index q in nonextensive statistical mechanics and 


with lim.y_,.i 1^].^ = (^). The Gauss binomial coeffi¬ 
cients can be arranged in a Pascal ^-triangle: 





























which reduces to the Pascal triangle for 7 —1. We define 
the Leibniz-like ^-triangle as: 


,(i) 

N^n,'y ' 


1 


1 




(57) 


with corresponding normalized triangle 


p(i) 

N,n,'y 


,(i) 

A^,n ,7 


fN\ ( 1 ) ’ 


and associated probabilities 


^( 1 ) 


(58) 


(59) 


In the limit 7 —>■ 1, triangles (57) and (58) reduce to the 
Leibniz triangle. 

Proceeding as in the previous Section, we can introduce 
the family 


„(i) 

(ly) _ Ar+2(i/'—1) ,71+1^—1,7 

^A/’,n,7 (1) ’ 


(60) 


deformation of the family ([^ , and the corresponding nor¬ 
malized version and its associated probability 




Indeed, let us start from 7 e (0,1): for 1, z/ > 1 

n 
N 


and ^ = V fixed. Then 


N,n^'y 


1 


1-7 


,(i) 

2 ( 1 /— 1,7 

n 

k^l 


^ _ ryN-{-‘ZV —1 


1-7^ 


X _ jk+N7]+v-l 

(i-7)nr=i 

^( 1 ) 

' 2(r/' —1),!/' —1,7 


0(1). (63) 


The relevant contribution in the probability distribution 
shape is therefore simply given by the binomial coefh- 
cient, and therefore, for N ^ 1, we recover the Gaussian 


distribution (see Fig. 4a) for all values of 7 G (0,1): 


4u) Af»l 
PN,un -- 


-2Ar(r,-i) = 

ttN 


(64) 


For 7 > 1 we can evaluate the asymptotic distribution 
as well. We can write 


1 


~(i/) 

PN:n,j 


1 ^ 

^^/ifc(n,7V;7), 

\n/ k—0 


(65) 


where 


fN\ 


y-fK—n L—'y 

ykJilj^l l-yN-k + i'-l+j 

hk{n,N-'y) ioi k = 

fN\ -rrn-fc 


for k > n, 
n, 

for k<n. 

( 66 ) 

For iV 1, ^ = 77 , and denoting hy ^ = k, we have 
that 


.(")n;r 


hNK{Nr],N;y) ~ 




(67) 


It is easily seen that the quantity J2nk finite 

only for 77 = 0 or 77 = 1, otherwise at least one of its 
addends diverges and ^ 0, n= I,... ,N — 1. 



'^N,n,y 

(61a) 


' A/’,n,7 

•sr^N N\ ’ 

Z^fc=0 \k)'’^N,k,-i 

»(i/) 

Pn,u. 

"( 33 ) 

Pn,u,j 

- \ri) 

(61b) 

V. 


Finally we get (see Fig. 4b for a numerical comparison) 

— -^-• ( 68 ) 


Remarkably, the y-triangles (61a) are neither strictly nor 


asymptotically scale-invariant since 


lim 


'^Ar-l,n.7 


A7—>-oo 

-^=77 fixed N,n,'y ' N,n-\-l,'y 


3 f”7<l, 

0 for 7 > 1 . 


In addition, probabilities (59) do not approach q- 


Gaussians with q 1 ^ as limiting distributions for large 
values of N. 


V. FINAL REMARKS 

We introduced three deformations — the a-, the /3- 
and the 7 -triangles — of the family of triangles pro¬ 
posed in [ 20 ], in order to analyze the robustness of the 
q-Gaussian family as attractors. Each one of the three 
proposed deformations depends on a single parameter in 
such a way that the undeformed family is recovered when 
the value of that parameter equals 1. We observed that, 
in all considered cases, the limiting distribution changes 
abruptly with respect to the undeformed case when a 
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(a) Probability distribution ^ for u = 2 and (b) Probability distribution 3 for different 

different values of 7 and theoretical prediction for values of ^ 

N = 10“^: not all data are represented for sake of 
clarity. 


FIG. 4: Numerical results for 7 -triangles. 




(a) Function qa{2) and its discontinuity for a = 1, (b) Function 5 ^ 3 ( 2 ) and its discontinuity for /3 = 1, 

see Eq. | |29[ |. see Eq. l |48[ l. 

FIG. 5: Functions qaii') and qp{v) evaluated for = 2 and different values of a and (3 respectively. 


deformed triangle is considered. Moreover, asymptot¬ 
ically equivalent deformations of natural numbers (the 
ct-deformation and the /3-deformation) exhibit different 
limiting behaviors, thus illustrating the high sensitivity 
of the limiting distribution with respect to the exact form 
of the deformation. However, for the considered asymp¬ 
totically scale-invariant deformations, the limiting dis¬ 
tribution of the deformed triangles is still a q-Gaussian, 
with a different value of q generically depending on the 
parameter of the deformation. Moreover, a discontinuity 
appears in q as function of the parameters a, /3 of the 
deformations, in correspondence of the undeformed case 
a = P = 1 (see Fig. [^. In particular 


lim -y— 

1 - qa[v) 


lim 


/3^i 1 - qfi{v) 


1 

1 - qi{v) 


-bl. (69) 


Remarkably, this discontinuity appears when we switch 
from a scale-invariant triangle to an asymptotically 
scale-invariant triangle. For /3-triangles, and similarly 
for a-triangles, this discontinuity expresses the fact that 

To exemplify this, let us introduce the following function: 






N 




G) 

■Pn,, 


n—0 


(71) 


The quantity is defined in such a way that, for 

iV — 00 , it remains finite. Indeed, Np^l^\p = Oil) in 
the -)■ 00 limit, and therefore I]n=o IPtv.u,/? = 
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FIG. 6 : Values of A^^(/3) for i/ = 2 and different values 
of N: the limiting function lim^_>oo A^^(/3) is also 
represented, including its discontinuity for /3 = 1. 


0(A). In Fig. @ it can be seen that the convergence 

of a|^^(/?) to A^c^\l3) := fimAr_>.oo A^^(/3) is not uniform 
and that a discontinuity for /3 = I appears in the V —> oo 
limit. 


The Qf-triangle is strictly stable under the action of the 
a-deformation for a ^ 1 , whereas a dependence of the 
limiting value of g on /3 (for /3 < 1) appears when the 
/ 1 -deformation is considered, although the limiting dis¬ 
tribution is still a g-Gaussian (Marsh et al. |T7] analyzed 
a probabilistic model in which gent bas a similar behav¬ 
ior with respect to a certain parameter of the model). 
The structure of relations (31) and (52) is quite com¬ 
mon in the literature: it has been observed, indeed, that 
for many systems characterized by a set of values of g, 
{gm}m 6 M, M := {0,1,... } C Z, a permutation of the 
indices m can be found such that Vm S M and a G (0, 2] 


a 

1 gm 


a 

- - 1 - m. 

1 - go 


(72) 


Finally, using the y-deformation, that is not an asymp¬ 
totically scale-invariant deformation, for 7 > 1 we ob¬ 
tain a limiting distribution that is not a g-Gaussian dis¬ 
tribution, and for 7 < 1 we always obtain a Gaussian. 
This fact suggests that the (asymptotic) scale-invariance 
property plays a central role in the robustness of the 
set of g-Gaussian distributions as limiting distributions. 
More specifically, the set of g-Gaussians appears to be ro¬ 
bust under asymptotically scale-invariant deformations. 
It may be interesting to investigate further the role of 
asymptotically scale-invariant deformations in the sta¬ 
bility of the g-Gaussian limiting distributions, to prop¬ 
erly identify the conditions under which the basin of g- 
Gaussians is an attractor for these probabilistic models. 
This fine of research could ultimately provide a deeper 
understanding of why there are so many g-Gaussians and 
g-exponentials in Nature. 


Appendix A: Asymptotic behavior of a- triangles 


In this Appendix we prove Theorem on the 

asymptotic distribution of a-triangles. We restate the 
theorem here for the reader’s convenience. 


Theorem. The family of triangles (20 1 with v G 
(1, -l-oo) and a G \ {0} satisfies the property 


N 


da,l 


4 v ) p ('^'1 

TN,n,a 


(Al) 


where we have introduced the properly centered and 
rescaled variable 


X := 






In Eq. (|A1[), Pq^(v){x) is a q-Gaussian with 


(A2) 


qa{iy) := 1 - 


u — 5, 


‘a,l 


1 - p fora^l, 
l-ph[ fora=l. 


(A3) 


Proof. The strategy of the first part of the proof consists 
in a generalization of the argument in |^. We want to 
evaluate the asymptotic behavior, for V ^ 1 , of 



(G 

'^Nn 


(A4) 


with rj := ^ fixed. Let us firstly observe that 


7 ,) _ r{2i^)r{N + I^) _ T{2v) 


E 


(1 - v)}. 


r{i^)T{N+ 2iy) ^ kl{N+ v + k)’ 

(A5) 


where we used the formula for the beta functior 0 


B(a, b) 


r(«)r(b) 

r(a -I- b) 


1 


0 


p-\l-tf-^dt = Y^ 

k=0 


{l-b)k 
k\{a k) ’ 

(A 6 ) 


This formula can be obtained using the Newton’s generalized bi¬ 
nomial theorem on the expression B(a, 6) := — 
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and {a)k ■= a(a + 1 )... (a + k — 1 ) is the rising factorial. 
The previous formula allows us to write down a general 
expression for the element as follows: 


'-Nn 




Ei-i)' 

2^0 
n oo 

EE 

i—O k—0 


AT—2,0 


k\ 


1 


ijN + v — i + k 

(A7) 


To evaluate the large N behavior, we use the saddle point 
approximation, 




i =0 


1 


ijN + v — i + k 


OO 

0 

where r] := ^. Inserting the previous term in the com¬ 
plete expression we have 

~ ^ (A9) 

Moreover, using the Stirling formula, 

^v(")r<;,l^h|l|(l-9).,r‘. (AlO) 

To complete the proof, we only need to evaluate the 
asymptotic behavior of „ for large N at fixed rj. 

Obviously, limQ,_j.i a ~ a 7 ^ 1 we have that, 

expression 

1V77(1 - 77) 


by direct computation on the expression (221 

(i^) iV»l 




(All) 


up to a multiplicative constant depending on a and 
ly. The thesis follows straightforwardly after a proper 
change of variable, ri 1 — ^ - 1 - and normaliza- 

tion. Observe also that the thesis holds for all real values 
1/ > 1. □ 


Appendix B: Asymptotic behavior of /3~triangles 


In this Appendix we prove Theorem in on the 
asymptotic distribution of /3-triangles. We restate the 
theorem here for the reader’s convenience. 


Theorem. The family of triangles (43b I with 1 / € IN and 
/3 > 0 satisfies the property 


N 


(u) p (p 

PN,n,p ^qi3(.ty)\-^) 


(Bl) 


2\/i^-x(/3) 
where we have introduced the function 

r .A [0 for (d>l, 

x(/3) := l-bd; 3 p-max J 1 , - I = J 1 for 13 = 1, 

^ [l-i /or 0</3<l, 

(B2) 

and the properly centered and rescaled variable 

(B3) 




In Eq. (Bl), is a q-Gaussian with 


qp{v) = 1 - 


' - x(^) 


1 - 

1 

u 

for P > 1 , 

1 - 

1 

u -1 

for (3=1, 

1 - 

0 

for 0 < /3 < 1 

Pv+l -0 


(B4) 


Proof. The proof of the Theorem strictly follows the 
proof of Theorem III.l in the only difference being 
the evaluation of the asymptotic behavior of the p 
coefficient for large N at fixed 

For /3 = 1, we have simply = 1 so there is 

nothing to do. 

For /3 7 ^ 1 and v = 2,3,..., denoting N — n =: Nrj, we 
can write 


!/-l 

k-N^,n,/3 — n ^k+v{P)~ 




TNr]—l-\-v ■ 


k=0 nfe=0 f^k+iJ+N{l-rj) 

Observe that the first product 

i^-i 




(B5) 


Bo bk+iy{/3) 


(B 6 ) 


fc =0 


is only a global factor not depending on N, n. We need 
to perform our asymptotic analysis only on the fraction 


UkS'h+.W) 


Uk^O bk+.+Nil-r,m 

We distinguish the two cases, 0 < /3 < 1 and /3 > 1. 
For 0 < /3 < 1 

Mvl,5 [77(1 - 77)]? , 


(B7) 


(B 8 ) 


where B_ is a certain constant. Indeed we have that the 
asymptotic behavior of the denominator of the fraction 
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(B7| can be evaluated as 

1 

Nrj-l + v , 

n i-s 


Eq. (B8) follows directly identifying i?_ = Bi[( 3 )B 2 {( 5 ). 


For /? > 1 we have instead that 


1 


pNr]-\-v 




(3 k + v + A^(l — rf) 


B+Nri{l - r]), 


(B15) 


+ pj Y (N (I — 1 ]) + ly) {l — r])P 
pN7^+‘'r{N + 2v) Y -T]) + i^- 


where -B+ is a certain constant depending on f3. Indeed, 
in Eq. (B5) we can write 


j^N-q+v 

(B9) 




Moreover, we have also 

7V?7—1 

n bk+u{p)= 

k^Q 

= ^ ^k+vjP) ^ JJ- bk+vjP) 


A^77 —1 Nr]—l+v 

^\{bk+AP) n (i- 

A^?7—1 


A;=0 

N{1 — rf) + V — 1 


k + V + Nil — rf) 


N + 2v-I 


Nr)—1 

bk+u{l3) ~ (I- 77 ) n 


- Li 


In the previous expression we introduced 

_Jl for 1 < fc < iV??-I, 

k+AP) ■ j ^ 

In the large N limit we have that 

Nq-l 00 , 

n 


. (BIO) 


(BlI) 


k=0 


k=0 


(B16) 


The remaining product can be written similarly as 


Nq-l Iwq-L , V 

n bk+A^) ^ B^{f3) n h+ j^Bi(/3)B2(I)iV77, 

fc^O k^O ^ ^ 

(BIT) 

where we used the definitions ( |B14 ). Eq. (B15) follows 
directly imposing B+ = Bi{f3)BAl)- 

Summarizing, for /3 > 0, ^ I, up to a global multi¬ 

plicative constant, 

« ~ [^^(1 - . (bI8) 


A^?7—1 


/c=0 


iV>l 




Nr}—1 


l\V' P{k + A 

r(i/ + i) ^{Nv + A 
^ BA/3)B2{l3)/3^m^r^K (B12) 

where we have introduced 

1 for I < fc < Nr] — I, 


Defining the function x(/?) as in (B2), we have 


~ [^(1 - ■ (B19) 


h+AP) ■ I ^ ^ ^ 

and the constant^ 


BA/3) := n 

k^O 

B2{/3) := 


1 _L _ 1-/3 

' fc+iv- 1-/3^+*^ 

min{/3,1} -h ^ ’ 


r(zz) 


'(" + ?)''' 


(B13) 

(B14a) 

(BI4b) 


Introducing now the variable 


:= 2^/1/-x(/ 3 ) (^1?- 


(B20) 


and properly normalizing we obtain the thesis. □ 
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